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Abstract 

Double-tag events in two-photon collisions are studied using the L3 detector at 
LEP centre-of-mass energies from ^/s = 189 GeV to 209 GeV. The cross sections 
of the e + e~ — > e + e~hadrons and 7*7* — > hadrons processes are measured as a 
function of the product of the photon virtualities, Q 2 = \JQ\Q\-, of the two-photon 
mass, W 77 , and of the variable Y = ln(W 2 /Q 2 ). The average photon virtuality is 
(Qi) = (Ql) = 16 GeV 2 . The results are in agreement with next-to-leading order 
calculations for the process 7*7* — > qq in the interval 2 < Y < 5. An excess is 
observed in the interval 5 < Y < 7, corresponding to W 77 greater than 40 GeV. 
This may be interpreted as a contribution of resolved photon QCD processes or the 
onset of BFKL phenomena. 



Submitted to Phys. Lett. B 



1 Introduction 



This letter presents a measurement of cross sections of two-photon collisions: e + e~ — > e + e~hadrons 
obtained with the L3 detector [I] using double-tag events, where both scattered electrons are 
detected in the small angle electromagnetic calorimeters. The virtualities of the two photons 
are defined as Q\ = 2EiEb(l — cos 6^), where Ef, is the beam energy, Ei and 6i are the mea- 
sured energy and scattering angle of the detected electron (i = 1) or positron (i = 2). The 
centre-of-mass energy of the two virtual photons, W 77 , is related to the e + e~ centre-of-mass 
energy, s/s, by Wj? ~ syiy 2 , with y i = 1— (Ei/E b ) cos 2 (#j/2). This is a good approximation 
in the kinematic range covered by this study, where W 2 7 is usually much larger than Qf. It is 
convenient to define the dimensionless variable Y: 

F = ln ^?> Q 2 = \foiQl (i) 

which depends mainly on the angles of the scattered electrons and allows the combination of 
the data collected at different values of y/s. 

Taking advantage of the good energy resolution of the small angle electromagnetic calorime- 
ters, W 77 is calculated as the missing mass of the two scattered electrons, W ee . This avoids an 
unfolding procedure to calculate W 77 from the effective mass of the hadrons observed in the 
detector, W V i S , which is the dominant source of systematic uncertainty on the measurement of 
the e + e~ — ► e + e~ hadrons cross sections for untagged 0E! and single-tag |HEj events. However 
the variable is more strongly affected by QED radiative corrections than W V i S . 

In perturbation theory, the cross section of the 7*7* — ► hadrons process is described in terms 
of a fixed order expansion in the strong coupling constant, complemented with the DGLAP [6 
evolution of the parton density of the photon. All two-to-two leading order (LO) processes, such 
as 77 — > qq (QPM) or, for example, jg — > qq or 7q — > gq (single resolved) and gg — > qq (double 
resolved), are taken into account in the Monte Carlo generators used to analyse the data. If 
the virtualities of the two photons are large and comparable, LO processes are expected to be 
suppressed relative to diagrams where multiple gluons are exchanged between the qq dipoles [7] 
coupling to each virtual photon. Examples of possible diagrams are given in Figure In 
leading logarithmic approximation, the resummed series of perturbative gluonic ladders can be 
described by the BFKL equation [B], which predicts a rise in cross sections as a power of W 77 , as 
if a "hard Pomeron" jU] was exchanged. The cross section measurement of two virtual photons 
is considered as a "golden" process to test BFKL dynamics JU]- After our first publication 
on the double-tag data at y/s = 91 GeV and 183 GeV [TT], an effort was made to improve 
the QPM calculation by including QCD corrections ^2]- The effects of varying the charm 
mass and the strong coupling constant were studied as well as the contribution of longitudinal 
photon polarization states [13J. The relative importance of perturbative and non perturbative 
QCD effects was also addressed by considering Reggeon and Pomeron contributions [T^ ITo']. 
There are also many efforts to include next-to-leading order (NLO) corrections in the BFKL 
model [IE]. 

The data, discussed in this letter, were collected at y/s = 189—209 GeV and correspond to an 
integrated luminosity of 617 pb _1 , for a luminosity weighted centre-of-mass energy 197.9 GeV. 
The observed value of Q\ is in the range 4 — 44 GeV 2 with an average value of (Qj) = 16 GeV 2 . 
The kinematic region E\^ > 40 GeV, 30 mrad < # 12 < 66 mrad and H^ 77 > 5 GeV is investi- 
gated. A study of asymmetric double-tag events (Q\ 3> Q\) at y/s = 91 GeV was previously 

^ Electron stands for electron or positron throughout this paper. 



2 



reported |T7j . 



2 Event Generators 

Two Monte Carlo generators, PHOJET [TH] and TWOGAM [T^j, are used to simulate double- 
tag two-photon events. Both use the GRV-LO [20] parton density in the photon and include 
all two-to-two LO 77 diagrams. They describe well single-tag events [3]. 

PHOJET is an event generator for pp, 7p and two-photon interactions, based on the Dual 
Parton Model. The electron-photon vertex for transversely polarized photons [21] is simulated. 
A transverse momentum cutoff of 2.5 GeV on the outgoing partons is applied to separate soft 
from hard processes [22] • PHOJET gives also a good description of untagged 77 — > hadrons 
events j2]. The electromagnetic coupling constant, a em , in PHOJET is fixed to the value for 
on-shell photons. 

TWOGAM generates three different processes separately: QPM, QCD resolved photon pro- 
cesses and non perturbative soft processes described by the Vector Dominance Model (VDM). 
The normalization of the QPM process is determined by the quark masses (m u = = 0.3 GeV, 
m s = 0.5 GeV and m c = 1.6 GeV), that of the VDM process is fixed by our measurement of the 
cross section of real photons [2|, while the normalization of the QCD contribution is adjusted 
to reproduce the observed number of data events. TWOGAM was recently upgraded to take 
into account QED soft and hard radiation from initial (ISR) and final state (FSR) electrons. 
The accuracy of the implementation of QED radiative corrections is checked with the program 
RADCOR [2S], using the channel e + e~ — > e + e~ . 

The data are mainly sensitive to initial state radiation which modifies the shape of the Y 
spectrum. Since the various processes have different Y dependences, the radiative correction 
affects them differently, as shown in Figure [2k- Here the cross sections are calculated in the 
generator level within the kinematic region defined above. The variables Q 2 and are 
calculated from the kinematics of scattered electrons. The relative contributions of QPM, 
VDM and QCD, as predicted by the TWOGAM program, including QED radiative effects, are 
given in Figure [2b and listed in Table ^ The VDM contribution is small and almost constant 
in our kinematical region. The resolved photon QCD contribution is negligible at low values of 
Y and increases to about 50% at high values. 

The dominant backgrounds are e + e~ — > e + e _ r + r~ events, simulated by JAMVG [23], an d 
single-tag two-photon hadronic events, where a hadron mimics a scattered electron. Other back- 
ground processes are simulated by PYTHIA 23] (e + e~ — > hadrons), KORALZ j2H] (e + e~ — > 
t+t-) and KORALW [23 (e+e~ -> W+W"). 

All Monte Carlo events are passed through a full detector simulation of the L3 detector 
which uses the GEANT j2B] and the GEISHA [22] packages and are reconstructed in the same 
way as the data. Time dependent detector inefficiencies, as monitored during the data taking 
period, are also simulated. The effect of the detector on the generated value of Y, Y gen , is 
presented in Figure [2t, where the distribution of value reconstructed from the hadronic system, 
Y vis , is shown in comparison with the quantity Y ee obtained from scattered electrons. The 
distortion and limited range of the Y vis spectrum, due to the effect of undetected particles, is 
evident. 
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3 Event Selection 



Double-tag two-photon events are recorded by two independent triggers: the central track 
trigger [3U| and the single- and double-tag energy triggers jHI] leading to a total trigger efficiency 
greater than 99%. 

Two-photon hadronic event candidates, e + e~ — ► e + e~hadrons, are selected using the fol- 
lowing criteria: 

• There must be two identified electrons, forward and backward, in the small angle elec- 
tromagnetic calorimeters. Each electron is identified as the highest energy cluster in one 
of the calorimeters, with energy greater than 40 GeV. The scattering angles of the two 
tagged electrons have to be in the range 30 mrad < < 66 mrad. The opening angle 
between the momentum vectors of the scattered electrons must be smaller than 179.5°, to 
reject Bhabha events. Figure El shows the distributions of Ei/Eb, Q 2 , 9i and \og(Ql/Ql) 
for scattered electrons. TWOGAM describes the shape of the distributions of 8i and Q 2 
better than PHOJET. 

• The number of particles, defined as tracks and isolated calorimeter clusters in the polar 
angle region 20° < 9 < 160°, must be greater than 5. The tracks are selected by requiring 
a transverse momentum greater than 100 MeV and a distance of closest approach, in the 
transverse plane, to the interaction vertex smaller than 10 mm. Isolated energy clusters 
are required to have energy greater than 100 MeV and no nearby charged track inside a 
cone of 35 mrad half-opening angle. 

• The visible hadronic mass W V i s , calculated from the four- vectors of all measured par- 
ticles, must be greater than 2.5 GeV in order to exclude beam-gas and off-momentum 
electron backgrounds. The distributions of W V i S and of the corresponding variable Y„i S = 
ln(W^ s / a/Q 2 ) are compared to Monte Carlo distributions in Figure and b. 

After these requirements, 491 events are selected with an estimated background of 49 
misidentified single-tag events and 32 events from the process e + e~ — > e + e~r + r~. Other 
background processes are estimated to contribute 6 events. The variable W 77 and the corre- 
sponding value of Y are calculated from the scattered electron variables, W ee and Y ee , shown 
in Figure and d. Good agreement is observed with both Monte Carlo generators. 

4 Results 

The differential cross sections of the e + e~ — > e + e~hadrons process with respect to the variables 
Q 2 , W 77 and Y are measured in the kinematic region: 

• E lj2 > 40 GeV and 30 mrad < 6» lj2 < 66 mrad 

• W 77 > 5 GeV 

The ranges 10 GeV 2 < Q 2 < 32 GeV 2 , 5 GeV < W 77 < 100 GeV and 2 < Y < 7 are indepen- 
dently investigated. The cross sections are derived in each interval as: 



where AN is the background subtracted number of events, £ is the total integrated luminosity 
and e is the selection efficiency. This is the ratio of the selected number of Monte Carlo events 
after the full detector simulation to the generated number of Monte Carlo events, including QED 
radiative corrections. An additional multiplicative factor £, discussed above and presented in 
Figure Et ; corrects the effect of QED radiative corrections. The results with and without this 
correction are given in Table E] for different bins together with the number of observed events 
and the selection efficiencies. The size of QED radiative corrections is estimated by TWOGAM 
using the relative proportions of the three components after adjusting the QCD component to 
the data. 

The systematic uncertainty on the cross sections due to the selection is 5%. It is dominated 
by the effect of a variation of the multiplicity cut from 4 to 6 particles. The uncertainty 
from the background estimation of single-tag events is 3.5% and that due to Monte Carlo 
statistics amounts to 1%. The uncertainty due to Monte Carlo modelling is estimated as 
6.4% by comparing PHOJET and TWOGAM without QED radiative corrections. To check 
the implementation of QED radiative corrections, the TWOGAM predictions for the e + e~ — > 
e + e~/i + yU _ process are compared to those of RADCOR. The difference is within 3% which is 
included as a systematic uncertainty. The different systematic uncertainties are summarised in 
Table El The different contribution from QPM, VDM and QCD as function of Y and gives 
an additional systematic uncertainty. A 20% variation of the QCD component results into an 
uncertainty of 0.3% at low values of Y and VF 77 and of 5.7% at large values. This uncertainty 
is about 0.5% over the full Q 2 region. 

The e + e~ — > e + e~hadrons cross sections after the application of QED radiative corrections 
are compared in Figure El to the PHOJET Monte Carlo and to LO and NLO calculations of 
7*7* — > qq [12J. In these calculations the mass of quarks is set to zero and a em is fixed to the 
value for on-shell photons. The predictions of these models are also listed in Table EJ These 
calculations describe well the Q 2 dependence of the data. For the W 77 and Y distributions, the 
QPM calculations describe the data except in the last bin, where the experimental cross sec- 
tion exceeds the predictions. Such an excess is expected if the resolved photon QCD processes 
become important at large Y, as illustrated in Figure Eb- The predictions of PHOJET, which 
includes the QPM and QCD processes in the framework of the DGLAP equation, also describe 
the data. A similar behaviour may also be obtained by considering the "hard Pomeron" con- 
tribution in the framework of BFKL 15J theories, while LO BFKL calculations were found to 
exceed the experimental values by a large factor jllj . 

From the measurement of the e + e~ — > e + e~hadrons cross section, a ee , we extract the two- 
photon cross section, a 7 » 7 *, by using the transverse photon luminosity function, L TT [2TJ 
E2L fee = L TT <7^*^*. a 7 * 7 » represents an effective cross section containing contributions from 
transverse (T) and longitudinal (L) photon polarisations: 



y*7* = CT TT + €\<T LT + t 2 UTL + Cl^LL + ^ClC2 J T TT COS 2(f<l{p - 4^i^ 2 J T TS COS <pd0 



(3) 



with 



Ci ~ Vi ~ e i = T— ^ ~~^~v? ' wllen Vi ^ 1 ( 4 ) 



i + (i-y-. 



2 

1 ) 



where (p is the angle between the e + e~ scattering planes in the two-photon centre-of-mass 
system. Using the GALUGA Monte Carlo program [22], the contribution of the interference 
terms ttt an d tts is found to be negligible for the QPM contribution, when Y > 3. In the 
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kinematical region studied, the average value of q is about 0.95. The experimental values of 
(j.y*-y* are presented in Table El and Figure |B1 in the same ranges considered above with and 
without QED radiative corrections. The measurements as a function of Q 2 are fitted by the 
form / = A/Q 2 , expected by perturbative QCD [rUl ll4j. The fit reproduces the data well with 
A = 81.8 ± 6.4 nb/GeV 2 and x 2 jd.o.f = 1.2/3. The average value of <t 7 * 7 . in the kinematical 
region considered is 4.7 ± 0.4 nb. The NLO calculations predict a decrease of (jry*™* as a 
function of or Y, which is inconsistent with the measurements at large values of and 
Y. 
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Table 1: Fractional contributions of the three processes, QPM, VDM, and QCD in different 
Q 2 , W J7 and Y intervals as predicted by the TWOGAM Monte Carlo including QED radiative 
corrections. 



11 



AD 2 


l(l 2 \ 






JDeiUie XVdQlclLl Vc CUIX eC IXUIlo 


ixnei xvaciiaLive cunecixuiio 


l uev ) 






6 


Arr lf\D 2 (r\h>ICp\I 2 \ 




10 — 14 


12.0 


128 5±12 4 


0.58 


0898 ± 0087 ± 0081 

• KJ KJ -J KJ — 1 — KJ , KJ KJ KJ 1 — 1 — UiUUUl 


0718 ± 0070 ± 0061 ± 0022 


14 — 18 


15.9 


102 0±11 2 


0.68 


0612 ± 0067 ± 0055 

KJ • KJ KJ -L j-J — 1 — \_l •\_I\J\J 1 — 1 — KJ • KJ KJ KJ KJ 


0522 ± 0057 ± 0044 ± 0016 


18 - 24 


20.5 


81 3± 9 8 


0.74 


0298 ± 0036 ± 0027 

• KJ £-1 -J KJ — 1 — KJ • KJ KJ KJ KJ — 1 — KJ • KJ KJ £-1 1 


0273 ± 0033 ± 0023 ± 0008 


24 - 32 


27.0 


24 8± 5 5 

^^T.t\J 1 J ' J 


0.77 


0065 ± 0014 ± 0006 

u. uuvu — i — vx • v/ vx j~ ~ — i — kjiKjkjkjkj 


0066 ± 0014 ± 0006 ± 0002 

KJ • KJ KJ KJKJ — 1 — KJ • KJ KJ 1~ a — 1 — KJ • KJ KJ KJ KJ — 1 — KJ • KJ KJ KJ -LJ 


AM/ 

i_i >yy 


IW \ 
\ vv ~i~i) 






JDclUIt! XVdjLlldjLl Vc CUIIcCllUilo 


A ft flr l-< orli q t n to c r\ v var^i i rxxx c 
xT-llt;I IXdjCIlclLl Vt; CUIItJCLlUIlD 




l uc V 1 


J_; Vcllljo 


6 


Arr 1 AW (r\h /CloV") 


Arr /AW fnh/HpV, 


5 — 10 


7.2 


67.3± 8.7 


0.37 


0594 ± 0076 ± 0053 

KJ • KJ KJ \J ^ — 1 — Vx • KJ KJ 1 KJ — 1 — KJ • KJ KJ KJ KJ 


0747 ± 0096 ± 0063 ± 0023 

W>VJ 1 a 1 — 1 — KJ * KJ KJ is KJ — 1 — KJ • KJ KJ KJ U — 1 — KJ . KJ KJ £j*J 


10 - 20 


13.9 


135 4±12 6 

-L '.J " J . X 1 L w * 


0.66 


0332 ± 0031 ± 0030 

WiUUUil — 1 — KJ • KJKJKJ -L. — 1 — KJ iKJKJKJKJ 


0263 ± 0024 ± 0022 ± 0008 

KJ * KJ ^ KJ U — 1 — KJ * KJ VX ^ a — 1 — Vx • Vx Vx ^ — 1 — KJ • KJ KJ KJ KJ 


20 - 40 


27.9 


102 1±11 1 

X V / ^ > X 1 1 X > X 


0.72 


0114 ± 0012 ± 0010 

\j • v ' x x x i v / • v / v j x _ i v / > v ' v > j_ v / 


0062 ± 0007 ± 0005 ± 0003 

UtVJUUiJ — 1 — KJ m KJ KJ KJ 1 — 1 — KJ • KJ KJ KJ U — 1 — KJ • KJ KJ KJ KJ 


40 — 100 


61.6 


65 1± 9 8 

\j ' j * j — i — t/ • i—) 


0.67 


0026 ± 0004 ± 0002 


0014 ± 0002 ± 0001 ± 0001 

KJ * KJ KJ ± a — 1 — KJ * KJ KJ KJ ^ — 1 — Vx • Vx Vx Vx J- — 1 — KJ • KJ KJ KJ 1~ 


A V 


lv\ 






Before Radiative corrections 


After Radiative corrections 






Events 


e 


dfTee/dr (pb) 


dcTee/dF (pb) 


2.0-2.5 


2.2 


51.6± 7.9 


0.52 


0.322 ± 0.049 ± 0.029 


0.315 ±0.048 ±0.027 ±0.009 


2.5-3.5 


2.9 


115.6±11.4 


0.73 


0.258 ± 0.025 ± 0.023 
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Table 2: Number of events, selection efficiencies, e, and differential cross sections da(e + e _ — > 
e + e~hadrons)/dQ 2 , d<r(e + e~ — > e + e~ hadrons) /dW 77 and da(e + e _ — > e + e~hadrons)/dY. All 
measurements are given before and after applying QED radiative corrections. The first uncer- 
tainty is statistical and the second systematic. The third uncertainty represents the effect from 
QED radiative corrections, including the 3% from Table 3. 
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QED radiative correction 
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Table 3: Contributions to the total systematic uncertainties on the measured cross sections. 
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Table 4: Predictions of LO and NLO 7*7* — > qq calculations and the PHOJET Monte Carlo 
generator as a function of Q 2 , W 77 and Y . 
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Table 5: The two-photon cross section, <r 7 * 7 », before and after applying QED radiative cor- 
rections, as a function of Q 2 , W 77 and Y. The first uncertainty is statistical and the second 
systematic. The third uncertainty represents the effect from QED radiative corrections, includ- 
ing the 3% from Table 3. 
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Figure 1: Examples of diagrams contributing to the process 7*7* — > hadrons : a) QPM, b) and 
c) 0(a s ) QCD corrections to the QPM diagram, d) photon-gluon fusion, e) one-gluon exchange 
and f) multigluon ladder exchange. 
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Figure 2: a) QED radiative corrections as a function of the variable Y, for QPM, VDM and 
QCD processes separately; b) the relative contributions of QPM, VDM and QCD processes 
in the TWOGAM Monte Carlo with QED radiative corrections included and c) Y determined 
using W V i S or W ee compared to the generated value, Y gen . Lines in a) and b) are drawn to guide 
the eye. 
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Figure 3: Distributions of a) Ei/E b , b) Qj , c) 9i and d) \og(Ql/Ql) for scattered electrons. 
The data are compared to Monte Carlo predictions, normalised to the total number of events 
in the data. The background is mainly due to e + e~ — > e + e~r + T~ and misidentified single-tag 
two-photon hadronic events. 
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Figure 4: Distributions of a) the effective mass of the detected particles, W vis , b) Y vis , c) the 
missing mass of the two scattered electrons, W ee and d) the variable Y ee . The range of W V i s and 
Y vis is limited to low values due to particles which escape detection. The data are compared to 
Monte Carlo predictions, normalised to the number of data events. 
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Figure 5: The differential cross sections of the e + e~ — ► e + e~hadrons process, in the kinematical 
region defined in the text, after applying QED radiative corrections, as a function of a) Q 2 , b) 
W 77 and c) Y. The LO and NLO predictions [T2] for the process 7*7* — > qq are displayed as 
the dashed and solid lines respectively. The dotted line shows the prediction of the PHOJET 
Monte Carlo. 
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Figure 6: Cross sections of the 7*7* — > hadrons processes as a function of a) Q 2 , b) W 77 , and 
c) Y in the kinematical region defined in the text, after applying QED radiative corrections. 
The dashed line represents the fit to the data described in the text. The NLO predictions of 
Reference ^] for the process 7*7* — ► qq are displayed as a solid line. 



19 



